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ABSTRACT 

The origin of the chemical composition of the intracluster medium (ICM) is 
discussed in this paper. In particular, the contribution from Type la supernovae 
(SNe la) to the ICM enrichment is shown to exist by adopting the fitting 
formulas which have been used in the analysis of the solar system abundances. 
Our analysis means that we can use the frequency of SNe la relative to SNe II 
as the better measure than MF e ,SNia/^Fe,totai for estimating the contribution of 
SNe la. Moreover, the chemical compositions of ICMs are shown to be similar 
to that of the solar system abundances. We can also reproduce the sulfur/iron 
abundance ratio within a factor of 2, which means that the abundance problem 
of sulfur needs not to be emphasized too strongly. We need more precise 
observations to conclude whether ICMs really suffer the shortage problem of 
sulfur or not. 

Subject headings: intracluster medium (ICM): individual (Abell 496, Abell 1060, 
Abell 2199, AWM 7)— X-rays— Supernovae 
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1. Introduction 



Primary work on intracluster medium (ICM) enrichment has focused on iron, since 
this was the only element that was accurately measured in a large number of clusters, 
prior to the launch of Advanced Satellite for Cosmology and Astrophysics (ASCA, [Tanaka 



|et al. 1994p . However, ASCA X-ray observations of clusters of galaxies provide the first 



opportunity to determine the element abundances of O, Ne, Mg, Si, S, Ca, Ar, and Fe for 
the bulk of the intracluster medium (ICM) in clusters of galaxies QMushotzky 19951 ). This 



abundance ratio pattern in the ICM provided a unique tool to probe the origin of these 
heavy elements. In particular, a precise analysis of these elemental abundances was done for 
four brightest, moderate temperature clusters of galaxies (Abell 496, 1060, 2199, and AWM 
7). The mean abundances of O, Ne, Si, S and Fe were reported as 0.48, 0.62, 0.65, 0.25, and 
0.32, respectively, relative to the solar photospheric abundances QMushotzky et al. 1996|) . 



Loewenstein & Mushotzky (1996) tried to explain the observed abundance ratios of 
these elements with respect to Fe (since Fe is the most accurately measured element) by the 
calculated yields from Type II supernovae (SNe II). They concluded that the abundance 
ratio pattern of ICM is very similar to that of SNe II. However, they emphasized that 
a significant SNe la contribution to the Fe enrichment could not be ruled out. In fact, 
Ishimaru & Arimoto (1997) concluded that a SNe la iron contribution of 50% or higher to 
the ICM enrichment could not be ruled out and might be favored by the observations. In 
the end, Gibson et al. (1997) pointed out the large uncertainties of the adopted massive 
star physics and uncertainties in the resulting SNe II's chemical composition, which make 
it impossible to determine the contribution of SNe la to the iron enrichment of ICM. As a 
result, we have no way of estimating the contribution of SNe la to the ICMs. 

There is another problem pointed out by Mushotzky et al. (1996), Loewenstein & 
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Mushotzky (1996), and Gibson et al. (1997). It is the lack of sulfur in the ICMs. Gibson 
et al. (1997) insist that there would appear to be no escape from the fact that the sulfur 
abundances are at odds with a SNe II-dominated origin to the ICM iron abundance 
(although the uncertainties are large). This appears to be entirely independent of adopted 
SNe II yields flGibson et al. 1997|) . Quantitatively, the observed ratio S/Fe favors the ICM 



SNe la iron fractional contribution to be larger than 80%. On the other hand, O/Fe, Si/Fe, 
Mg/Fe, and Ne/Fe favor a ratio under 60%. 

In this paper, we adopt the fitting formula used in the analysis of the solar system 
abundances ( |l'sujimoto et al. 1995| ) and a \ 2 fitting formula to obtain the relative 



contributions of SNe la and SNe II to the enrichment of ICM. We anticipate that these 
formulas will become good measures for estimating the fractional contribution of SNe la. 
We will also adopt many models of SNe la (model W7 and WDD2, [Nomoto et al. 1997]) and 
SNe II ( Nagataki et al. 1997|) which take into account the effect of axisymmetric explosion, 



in order to reexamine the shortage problem of sulfur. 

In section ||, the results of X-ray observations of clusters of galaxies are presented. In 
section we show a range of calculated yields from SNe la and SNe II. In section |], we 
present our analysis and results. Summary and discussion are given in section |^. 



2. Observations 



We show the abundances derived from SIS (solid state imaging spectrometer) in 
Table [I]. [Mi/Fe] is defined as log(Mi/Fe)icM — log(Mi/Fe) so i a r- The error bars indicate 90% 
confidence intervals. We note that the meteoritic abundance scale of Anders & Grevesse 
(1989) is adopted. Oxygen, silicon, and iron have the best-determined abundances. On the 
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other hand, the accuracy for magnesium, argon, and calcium is very poor flMushotzky et aE 
19961) . 

EDITOR: PLACE TABLE | HERE. 



3. Models of SNe la and SNe II 



There is a range of SNe II models available in the literature. In the present situation, 
there are large differences in the predicted yields since there is a wide variety of input 
physics (e.g., criterion for convection, reaction rates, and the way to initiate the shock 
wave) ( [Aufderheide et al. 199 1| ; |Gibson et al. 1997|) . 



In order to compare the differences quantitatively, we make use of the formalism 
presented by Ishimaru & Arimoto (1997). The zth element yield averaged over the SNe II 
progenitor initial mass function (IMF), cf) oc m~ x , is written as 

_ Im" 2/i,SNii(m)0(m)m _1 dm 

< 2/i,SNii >- r m u 77 — \ — rTj • 

J m " (p{m)m Mm 

The bounds for SNe II progenitors are taken to be mi = 10M Q and m u = 5OM , 
respectively, and an IMF slope x = 1.35 QSalpeter 1955|) is adopted throughout. While this 



slope is of prime importance for arguments concerning the absolute mass of elements in the 
ICM (e.g., |Loewenstein fc Mushotzky 1996| ), it is less so for abundance ratios flGibson et al. 
P557D . 

In imitation of Gibson et al. (1997), we summarize < y^sNii > f° r a range of SNe II 
models in Table |^. We add the yields predicted by Nagataki et al. (1997) and explain it 
here. They performed 2-dimensional hydrodynamical calculations and studied the changes 
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of the chemical compositions using a large nuclear reaction network containing 242 nuclear 
species. Model N97S1 in Table |^ is a spherical explosion model, N97A1 and N97A3 are 
models of axisymmetric explosions in an increasing order of the degree of deviation from 
spherical symmetry. For the axisymmetric models N97A1 and N97A3, the initial velocity 
behind the shock wave is assumed to be radial and proportional to r x 1 + a ^> s ^ 29 ) ; w here r, 9, 
and a are the radius, the zenith angle and the model parameter that determines the degree 
of deviation from spherical symmetry, respectively. Nagataki et al. (1997) took a = for 
model N97S1, a = § for model N97A1, and a = J for model N97A3. The lar e:er a sets, 
the more asymmetric the explosion becomes. They used the same distribution also for the 
thermal energy. Half of the total energy appears as kinetic energy and the other half as 
thermal energy. 

As we can see from Table 0, almost all < y^sNii > of Nagataki et al. (1997) are in the 
range of uncertainties of theoretical predictions. However, the amount of sulfur tends to 
become smaller as a gets larger. In fact, N97A3 predicts the least amount of sulfur among 
all models. This tendency is expected to be good for the reproduction of the observed 
amount of sulfur in ICMs. 

We also show the SNe la yields < y^sNeia > at the lower portion of Table [|. W7 is the 
model of the simple deflagration and WDD2 is that of the delayed detonation ( Nomoto et 



al. 1997|) . As pointed by Gibson et al. (1997), all SNe la yields have been tied exclusively 



to the W7 model. In this paper, the uncertainty of SNe la yields is taken into consideration 
for the first time. 

EDITOR: PLACE TABLE | HERE. 



4. Analysis & Results 
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4.1. Previous Analysis 



At first, the analysis which is done in the previous papers is reexamined. We calculate 
the abundances ratios [M ; /Fe] for M ; = O, Ne, Si, and S as a function of M Fe ,sNia/^Fe,totai- 
^Fe,SNia/-^Fe,totai is the contribution of SNe la to the ICM iron enrichment and is given as 

M FeiSNIa _ r < l/Fc,SNIa > (0 < r < 1) (1) 

M F e >to tal T < ?/Fe,SNIa > +(1 - r) < ?/Fc,SNII >' ~ ~ 

where r indicates the relative occurrence frequency of SNe la. The result is shown in 
Figure [I]. Horizontal dot lines mean the average of ASCA SIS data. 

EDITOR: PLACE FIGURE Q HERE. 

[O/Fe] and [Ne/Fe] are not changed even when the model of SNe II is changed. 
This is because the abundances of O and Ne are not determined mainly by the explosive 
nucleosynthesis but by the nucleosynthesis during the stellar evolution. We also note that 
the abundance of Fe is mainly determined by the amount of 56 Ni, which decays to 56 Fe. 
Since the position of the mass cut is also determined by the amount of 56 Ni, Fe abundance 



is not changed among three models, too flNagataki et al. 19971) . On the other hand, [Si/Fe] 



and [S/Fe] are sensitive to the selection of the model of SNe II. This is because Si and S 
are mainly synthesized during the explosive nucleosynthesis. Si and S abundances tend to 
decrease along with the degree of deviation from spherical symmetry. Because of this, the 
[S/Fe] ratios give smaller M Fej s NIa /M FeiTotal when W7 and the axisymmetric models are 
adopted. 

As a result, the observations of O, Ne, Si, S, and Fe favor the ICM SNe la iron 
fractional contribution in the range 30% < M Fej sNia/^Fe,Totai < 60% for the models W7 
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and N97A3. On the other hand, it is in the range 30% < M Fe> sNia/^Fe,Totai < 90% for 
the models W7 and N97S1. Even worse, the averaged [S/Fe] ratio can not be reproduced 
when model WDD2 is adopted, although the theoretical [S/Fe] ratio with a proper ratio of 
^Fe,SNia/^Fe,Totai can be within the error bars. 

As we see, the ratio of M Fe gNIa /M FeiTota i is very hard to determine, which was also 
indicated in the previous papers ([Loewenstein fc Mushotzky 1996| ; [Ishimaru fc Arimoto 
19971 ; IGibson et al.~1997|) . 



4.2. g(() and \ 2 ~ fitting 



Next, we define ysm& as the sum of < y^sNia >, which is the heavy element mass of 
SNe la. ysNii is defined in the same way. We also define the abundance pattern Xj as 

Xi = C^,SNIa/?/SNIa + (1 ~ C)Z/i,SNIl/Z/SNII (0 < C < 1) ( 2 ) 

which is to be compared with the observed XjjcM- ^i.icm is defined as Zi/Y^iZi, where Zj 
is the observed abundance of the z-th element per unit mass. ( is the mass fraction of SNe 
la's matter in ICM. The relation between r and ( is written as 

l-T u u y u {l-0 U 
where and uu represent the mass fraction of heavy elements ejected into the interstellar 

gas from SNe la and SNe II, respectively. These values are estimated to be 0.27 and 0.22 in 

the solar neighborhood from the numerical calculation ( [Tsujimoto et al. 1995|) . 

The most probable value of ( = ( p is determined by minimizing the following function 
(Yanagida et al. 1990): 

n 

9(0 = X][log£i,iCM - \ogXi] 2 /n (4) 
i=i 
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where i runs over the heavy elements considered in the minimization procedure. We 
summarize those elements in Table [|. We will give another method to obtain £ p , that is, x 2 
fitting. As is well known, 

X 2 -J2( Xi ^~ x ^ 1CM \ 2 (5) 



obeys the x 2 distribution for n degrees of freedom. We note that Oi is the standard deviation 
for XijcM- Since we have no information about <7j, we use the error bars which indicate 90% 
confidence intervals ( [Ishimaru fc Arimoto 199~T| ) instead. 

We show the results for g(() and x 2 in Figure [|, £|, and ||. It should be noted 
that g(() and x 2 have a local minimum/maximum for each model, which means there is a 
favorable value of ( for the reproduction of the chemical composition of ICMs. Moreover, ( p 
is in the range of < ( p < 0.1 for almost all clusters, which represents that the contributions 
of SNe la to the ICM enrichment are similar among four ICMs. It should also be noted 
that ( p = 0.09 ± 0.01 is obtained for the solar system abundances QTsujimoto et al. 1995| ), 
which will mean that the chemical compositions of ICMs are similar to that of the solar 
system abundances. 

Now, there is one question: why can ( p be determined in the narrow range ( p < 0.1 in 
spite of the large uncertainty of the inferred ratio M Fe) sNia/-^Fe,totai ? The answer is shown 
in Figure || which gives the relation between £ and M Fe SNIa /M Fe total . As can be seen from 
the figure, M Fej sNiaA^Fe,totai is very sensitive to £ if £ is in the range of < ( p < 0.1. This 
is because the mass fraction of Fe in the SNe la's ejecta is quite high. As a result, it is 



impossible to determine M Fe ,sNia/^Fe,totai exactly, as we showed in section fO| . On the other 
hand, we can use ( as a better measure than M Fe SNIa /M Fc total to represent the contribution 
of SNe la to the ICM enrichment. 



EDITOR: PLACE FIGURE | HERE. 
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EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE § HERE. 

Finally, Figure [7] and [8] show the normalized abundance pattern {xi/x-p e )iCM/{xil '^Fe)sN 
with the most probable value ( = ( p which is obtained from the analysis of g(Q for each 
cluster. The abundance of sulfur seems to be reproduced better by N97A3 than N97S1, 
to be sure, however, the sulfur abundance can be reproduced within a factor of 2 even by 
N97S1. We note that the supernova abundances agree with the solar system abundance 
ratios within a factor of 2-3 for typical species ( [Hashimoto 1995| ), which will reflect the 
precision of the numerical calculations of supernova nucleosynthesis. Because of this, we 
think the sulfur abundance problem needs not to be emphasized too strongly. We wait for 
more precise observations to conclude whether ICMs really suffer the shortage problem of 
sulfur or not. 

EDITOR: PLACE FIGURE HERE. 
EDITOR: PLACE FIGURE E HERE. 



5. Summary and Discussion 



In this paper, the origin of the chemical composition in ICMs is discussed. At first, 
the analysis which has been done in the previous papers is reexamined. As expected, 
the ratio of M Fe) g NIa /M FeiTota i is very hard to determine. Quantitatively, the observations 
of O, Ne, Si, S, and Fe favor the ICM SNe la iron fractional contribution in the 
range 30% < M FCiS Nia/M FCiTota i < 60% for the models W7 and N97A3. The range is 
30% < M FCiSNIa /M FeiTota i < 90% for the models W7 and N97S1. Even worse, the averaged 
[S/Fe] ratio can not be reproduced when model WDD2 is adopted. Our analysis shows that 
^Fe,SNia/^Fe,Totai is not a good measure for estimating the contribution of SNe la at the 
moment to the ICMs. 

Next, the analysis of g(Q and x 2 fitting are done in order to investigate whether these 
can be better measures for SNe la's contribution. We have shown that g(() and x 2 nave 
a local minimum/maximum for each model, which means there is a favorable value of ( 
for the reproduction of the chemical composition of ICMs. Moreover, ( p is in the range 
of < ( p < 0.1 for almost all clusters, which means that the contributions of SNe la to 
the ICM enrichment are similar to each other and to the solar system abundances. This 
analysis indicates that ( is a better measure than M Fei sNia/-^Fe,totai- 

The reason why ( p can be determined in the narrow range ( p < 0.1 in spite of the large 
uncertainty of the inferred ratio M Fej sNia/^Fe,totai is that M Fei sNia/A^Fe,totai is very sensitive 
to C if C is in the range of [0, 0.1]. This is because the mass fraction of Fe in the SNe la's 
ejecta is quite high. 

Finally, we have investigated the sulfur abundance problem. The normalized abundance 
pattern (xj/x Fc )icM/(£«/£Fe)sN with the most probable value ( = ( p which is obtained from 
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the analysis of g(() shows that the sulfur abundance can be reproduced within a factor of 
2. Taking the uncertainty of the numerical calculation of supernova nucleosynthesis, we 
think that the sulfur abundance problem needs not to be emphasized too strongly. We need 
more precise observations to conclude whether ICMs really suffer the shortage problem of 
sulfur or not. The same applies for Ar and Ca. We hope that there will be a large number 
of more precise observations of ICMs to have a more advanced discussion. 

We would like to thank W. Hillebrandt for useful comments. We are also grateful to 
V. Saar for his kind review of the manuscript. This research has been supported in part 
by a Grant-in-Aid for the Center-of-Excellence (COE) Research (07CE2002) and for the 
Scientific Research Fund (05243103, 07640386, 3730) of the Ministry of Education, Science, 
and Culture in Japan and by Japan Society for the Promotion of Science Postdoctoral 
Fellowships for Research Abroad. 
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Fig. 1. — Theoretical ratio of relative abundances of heavy elements, [O/Fe], [Ne/Fe], [Si/Fe], 
and [S/Fe] as a function of the SN la fraction in the iron synthesis. Solid, short-dashed, and 
dot lines correspond to N97S1, N97A1, and N97A3 models, respectively. Left: W7 model is 
adopted for SNe la. Right: WDD2 model is adopted. 
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Fig. 2.— g(C) as a function of (. W7 and WDD2 are used for SNe la. Left: Abell 496. 
Right: Abell 1060. 
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Fig. 3.— Same as figure 2 but for Abell 2199 and AWM 7. 
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Fig. 4. — Confidence level (%)of x 2 fitting as a function of (. Left: Abell 496. Right: Abell 
1060. 
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Fig. 6. — Relation between ( and M Fe) sNeia/^Fc,Totai- Solid, dot, short-dashed, and long- 
dashed lines correspond to W7+N97S1, WDD2+N97S1, W7+N97A3, and WDD2+N97A3, 
respectively. 
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Fig. 7. — Normalized abundance pattern (xj/xF e )icM/( 3: i/2 ; Fe)sN with the most probable 
value ( = (p which is obtained from the analysis of g((). Left: Abell 496. Right: Abell 1060. 
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Fig. 8.— Same with figure 7 but for Abell 2199 and AWM 7. 
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Abell 496 Abell 1U6U Abcll 2199 AWM 7 

07Fe] +0.15(-0.15,+0.36) -0.08(-0.33, +0.08) -0.06(-0.45, +0.15) -0.01(-0.37, +0.20) 



Ne/Fe] +0.28(+0.10, +0.43) +0.15(+0.00, +0.25) +0.12(-0.12, +0.28) -0.16(-0.61, +0.07) 

Mg/Fe] -0.01 -0.43, +0.21) -0.31 -0.92,-0.06) -0.05 -0.48, +0.18) - (< -0.39) t 

Si/Fe] +0.06 -0.09, +0.18) +0.14 +0.07, +0.21) +0.26 +0.16, +0.33) +0.03(-0.08, +0.13) 

S/Fe] -0.23(-0.62,-0.02) -0.30(-0.55, -0.15) -0.31(-0.86, +0.28) -0.47(-0.97, -0.23) 

Art/Fe] -0.18(< +0.20) -1.31(< -0.23) - (< -0.31) - (< -0.52) 

Cat/Fe] -0.25(< +0.20) - (< -0.23) - (< -0.09) - (< -0.39) 

Ni/Fe] +0.49(+0.25,+0.65) +0.23(-0.33, +0.61) +0.37(+0.00, +0.54) +0.19(-0.19, +0.40) 



Table 1: Abundances of ICM normalized by the meteoritic abundances, f : Nuclei omitted 
in the analysis of g((). 



Yield Source 


< yo.sNii > 


< yNc.SNII > 


< 2/Mg,SNII > 


< ysi.sNii > 


< ys,SNii > 


< yFc,SNII > 


A96 


0.593 


0.101 


0.054 


n/a 


n/a 


0.071 


T95 


1.777 


0.232 


0.118 


0.133 


0.040 


0.121 


T95+M92 


0.923 


n/a 


n/a 


n/a 


n/a 


0.121 


W95;A;10- 4 Z Q 


0.806 


0.095 


0.036 


0.104 


0.059 


0.073 


W95;B;10q 4 
W95;A;Z 


1.455 


0.223 


0.066 


0.118 


0.065 


0.085 


1.217 


0.181 


0.065 


0.124 


0.058 


0.113 


W95;B;Z Q 


1.664 


0.265 


0.094 


0.143 


0.064 


0.141 


N97S1 


1.749 


0.229 


0.124 


0.139 


0.042 


0.086 


N97A1 


1.726 


0.231 


0.126 


0.106 


0.026 


0.086 


N97A3 


1.694 


0.233 


0.123 


0.098 


0.020 


0.088 



Yield Source 


< yO.SNIa > 


< yNc,SNIa > 


< 2/M K ,SNIa > 


< ySi.SNIa > 


< ?/S,SNIa > 


< yFc,SNIa > 


W7 

WDD2 


0.148 
0.069 


0.005 
0.0009 


0.009 
0.005 


0.158 
0.272 


0.086 
0.168 


0.744 
0.700 



Table 2: Average stellar yield in solar masses for the SNe grids. A96 = Arnett 1996, T95 
= Tsujimoto et al. 1995, M92 = Maeder 1992, W95 = Woosley & Weaver 1995, N97 = 
Nagataki et al. 1997c. N97S1, N97A1, and N97A3 stand for the model of the spherical 
explosion and those of the axisymmetric explosion, respectively. The degree of deviation 
from spherical symmetry is larger in N97A3 than N97A1. W7 is the model of the simple 
deflagration and WDD2 is that of the delayed detonation (Nomoto et al. (1997)). 



